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Abstract 

The procedure used for structure function measurements at HERA is briefly 
described and related to the properties of kinematic reconstruction. The recon- 
struction methods of the inclusive deep inelastic scattering variables are reviewed 
and their sensitivity to the energy and angle miscalibrations are discussed in de- 
tail. New prescriptions are introduced and related to the standard methods in 
order to optimize the F2 structure function measurement over the widest kine- 
matic range, both in the low x, low Q 2 and in the high x, high Q 2 regions. The 
prospects for the future high Q 2 studies are briefly discussed. 



1 Introduction 

The measurement of the structure functions of the nucleon is a major tool for the study 
of the strong interaction and the behaviour of the parton densities in the hadrons. To 
reveal possible non-standard small deviations from their well established behaviour 
described by the DGLAP [|J evolution equations, such as BFKL effects at low x or the 
presence of intrinsic charm in the proton at high x to name just two of them, requires a 
precise reconstruction of the deep inelastic scattering (DIS) kinematics over the widest 
possible kinematic range. With the advent of the HERA electron-proton collider, this 
reconstruction no longer needs to rely on the scattered lepton only, since the most 
important part of the hadronic system is visible in the almost hermetic HI and ZEUS 
detectors. This redundancy allows for an experimental control of the systematic errors 
and of the radiative corrections to the structure function measurement if it is based on 
several independent methods to determine the usual DIS kinematic variables x, y, Q 2 : 

x = ^- q y = £| Q 2 = -(k - k') 2 = -q 2 = xys (1) 

with s being the ep center of mass energy squared, P, k the 4-vectors of the incident 
proton and lepton, and k' of the scattered lepton. Since many different reconstruction 
methods have already been used at HERA, it is a natural objective to optimize this 
reconstruction, and to try to find "the best" method, or at least to justify the use of 
a given method instead of a kinematic fitting procedure for instance. In this report 
we briefly sketch in section 2 the procedure used so far to measure F 2 at HERA, and 
relate it to the effects of kinematic reconstruction. In section 3 we review the methods 
of kinematic reconstruction used at HERA, and we classify them on the basis of their 
properties. In section 4 we discuss possible improvements of these reconstruction meth- 
ods, in particular those related to the study of low x physics. In section 5 we study 
the effect of the hadronic final state and scattered electron reconstruction errors on the 
kinematic methods in order to understand the choices made in this field by the two 
HERA collaborations. In section 6, the high Q 2 case is treated in more detail due to its 
future importance and also since there are new possibilities in this kinematic regime. 
In conclusion, we briefly provide some prospects on the influence of these technical 
matters on the structure function program of the next decade. 



2 Structure Function Measurements 



The cross-section for the DIS reaction e + + p — > e + + X with unpolarized beams is: 
2T( " %F 2 (x,Q 2 )-y 2 F L (x,Q 2 )-Y_xF 3 (x,Q 2 )} (l + 6 r (x,Q 2 )) (2) 



dx dQ 2 xQ 4 

In this equation a is the electromagnetic coupling, F 2 is the generalized structure 
function which reflects both photon and Z° exchange, F L is the longitudinal structure 
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function, F3 is the parity violating term arising only from the Z° exchange, 5 r is the 
electroweak radiative correction. The helicity dependence of electroweak interactions, 
is contained in the functions Y±(y) = 1 ± (1 — y) 2 . 

At HERA, these three structure functions of the proton can be measured, in par- 
ticular F 2 which has already been measured over several orders of magnitude in x and 
Q 2 with a precision of about 5 %. Although apparently only one observable ( ^xIq 2 ) 
is related to the three structure functions, the problem can be solved because in some 
kinematic region only one of them has a relevant contribution and/or because the 
beam conditions (energy, lepton sign) can be changed, hence changing the coefficient 
in front of each of them, for a given d f d Q 2 ■ In the following, we will not study the 
procedure to derive them from the cross-section measurement (which are described in 
for instance), but concentrate on the determination of d f d q2 ■ 

Generally, the value of d f dQ 2 at the point (x a , Q 2 ) is experimentally determined 
by the number A^ ("D" stands for Data) of DIS events observed in an (x, Q 2 ) bin 
A centered around {x ,Q 2 ), normalized to the integrated luminosity L D accumulated 
during the data taking and corrected by an acceptance factor which depends on 
the event selection cuts and on the detector response: 

d 2 a C A Ng (3) 



dxdQ 2 L D 

Ca is a numerical factor which takes into account the surface of the bin. Actually, 
resolution effects cause migration of events from one bin to other bins, so the previous 
relation should be treated in fact as a matrix relation, with the cross-section being 
related to all bins via the inverse of a matrix T(ij), in which every matrix element (ij) 
gives the probability that an event originating from the bin i is reconstructed in a bin 
j. Thus the double differential cross-section measured at the center of the bin (z) is: 

(4 » 

The matrix T can be obtained by simulating a large sample of DIS events, provided the 
simulation is able to reproduce the detector effects. However to solve equation (4) one 
needs to invert this matrix. This is a non-trivial task due in particular to numerical 
instabilities j|, |, [|, which has not been done yet for the published HERA structure 
function measurements. Instead, an iterative procedure is used to solve equation (3). 
Obviously, if the structure function used in the simulation is equal to the one to be 
measured, and assuming that T£ = ("5" stands for simulation), we have: 



d 2 a A? L s d 2 a s 



A 



dxdQ 2 L D N s A dxdQ 2 



(5) 



This equation is solved iteratively, starting the simulation from a "guessed" structure 
function parametrization, and replacing it at each iteration by the parametrization 
obtained from a fit to the structure function "measured" at the previous iteration. 
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In three or four iterations the result is stable within one % or less (see for a 
more detailed discussion) in the measureable regions. These measureable regions are 
empirically characterized by values of Ta close to unity In order to define them in a 
more rigorous way and to understand their relation to the kinematic reconstruction, 
we will now study the Ta factors. 

We distinguish two ways of defining the bin A: one (A, t) based on the "true" 
kinematic variables, defined at the hadronic vertex; the other (A, r) based on the 
reconstructed variables. We also consider the effect of the "event selection cuts" on 
the number of events (N A ) compared to the number of events before the cuts (N&). 
These "cuts" are imposed to improve the precision of the measurement, but can have 
an influence on the distribution of the events in the kinematic plane, so the two cases 
must clearly be separated. With these definitions, we deduce from eq. (3) that Ta can 
be expressed as : 

Ta = -Tf-^ = tA,t ■ A A with e A ,t = -vr— and ^a = -^f- (6) 

The first term (eA,t) characterizes the "efficiency" of the cuts in the bin A, while 
the second one can be defined as the "smearing acceptance" of the bin A, since only 
the smearing of the kinematic reconstruction is involved in its variations. Since the 
cuts are choosen to have a high efficiency, we will not discuss here the difference at 
the percent level which might occur between €a as determined on the data and on 
the simulation, but focus on the behaviour of the smearing acceptance, which can 
have large variations (up to hundreds of percents) across the kinematic plane, hence 
determining the measureable regions. The smearing acceptance variations are better 
studied considering the (A, 1) subset which contains the events belonging to (A, i) and 
(A,r). Then A a can be rewritten as 

Aa = with Sa = -T^r 1 and P A = - ; N A ,i = N A ,tnA,r (7) 

-TA ^A,t ^A,r 

The Sa and Pa are refered as the stability and the purity of the bin A, since they 
characterize respectively: 

• the proportion of genuine events of a bin, which are reconstructed in the same 
bin. Sa characterizes the number of events which migrate outside of the bin A. 

• the ratio of genuine events of a bin reconstructed in the same bin divided by the 
total number of events reconstructed in that bin. This last number is influenced 
by the number of events originating from other bins which migrate into A. 

The obvious goal is to maintain a stability and a purity as close to unity as possible, 
hence A a will also be close to 1. However A a can be close to unity even for a low 
stability bin, if its purity roughly matches its stability. In order to ensure a reliable 
measurement of the cross-section, we will thus enforce conditions on Sa and Pa sep- 
arately, rather that only on their ratio. To make a concrete example let us examine 
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the three standard reconstruction methods used for DIS neutral currents at HERA, in 
the (x, Q 2 ) binning used by the HI collaboration, i.e 8 (5) bins per order of magnitude 
in Q 2 (x) ||. The smearing acceptance, stability and purity for the HI detector are 
shown in fig. [I] for the electron, Double- Angle |I| and S methods, in x bins, at Q 2 



20 GeV . As mentioned, the acceptance can be equal to unity, while the stability 



Smearing Acceptance 



0.5 



A e 
□ DA 
• I 



a n □ □ □ 



icr 3 

Q 2 = 20 GeV 2 



A • 

A □ • 
A • 

A A . 



10" 



Stability 



0.8 



0.6 



'0.4 



A e 
□ DA " 
• I . 



A 

A 
□ 



A • 

□ □ • « 

A 5 



A □ 
A „ 



10- 3 

Q 2 = 20 GeV 2 



Purity 



0.8 



0.6 



0.4 



0.2 



A A 
-A A 
. • • 



A e 
□ DA 
• I . 



Q = 20 GeV 



Figure 1: Smearing Acceptance, Stability and Purity at Q 2 =20 GeV 2 for the 3 kine- 
matic reconstruction methods: e (open triangles), DA (open squares), £ (closed circles). 

and the purity are low (~ 25%), like for instance at the lowest x point in the DA 
method. Conversely, by requiring a "reasonably" high stability and purity (we will 
use throughout this paper a minimum of 30%, as the HI collaboration did for this 
binning 0), we ensure an acceptance close to unity. 

We observe in fig. |l] that the stability has a more regular behaviour than the purity 
or the acceptance, since it reflects only the combination of the resolutions in x and Q 2 
in a given bin. The stability can be increased by enlarging the size of the bin. Although 
the stability and purity are related quantities (they share the same numerator), the 
purity is more irregular because it is influenced by events coming from other bins, which 
may be populated in a different way, and have different resolutions in the kinematic 
variables. Like the stability, the purity can also be increased by enlarging the bin size, 
however the migrations inside the bin depend on several factors (resolutions in different 
bins, population of different bins, structure function values in regions which are not 
"measurable", radiative effects), rendering the control of the purity more delicate than 
the stability. Thus, in the following, when comparing different methods, we will make 
the comparison on their purity, but only in the bins having a minimum stability (chosen 
to be 30%). 



3 Kinematic Reconstruction Methods at HERA 



In the naive quark-parton model (QPM), the lepton scatters elastically with a quark of 
the proton, and the two body final state is completely constrained using two variables, 
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if we know the initial energies labeled E D and P of the electron and proton. Similarly, 
the DIS variables can be determined using 2 independent variables, which can be 
the energy (E) of the scattered electron, its polar angle[](0) or independent quantities 
reconstructed out of the hadronic final state particles. For instance S, obtained as the 
sum of the scalar quantities E h —p zh of each particle (assumed to be massless) belonging 
to the hadronic final state, pr,h as its total transverse momentum or the inclusive angle 
7 of the hadronic systemF] which corresponds to the angle of the scattered quark in the 
QPM: 



s = J2( E h - Pz,h) Pr,h = J(^Px,h) 2 + (Y,Py,h) 2 tan- = (8) 



h 



2 PT,h 



Eh, Px,h, Py,h,Pz,h are the four-momentum vector components of each hadronic final state 
particle. £ is by construction minimally affected by the losses in the forward direction 
due to the beam pipe hole in which the target jet and the initial state gluon radiation 
tend to disappear, pT t h covers the other spatial dimensions, and is more sensitive to 
forward losses. Thus to improve the kinematic reconstruction we should avoid using 
Px,h directly. We can use 7 instead, which carries the px,h information and is better 
measured since in the ratio E/p^/j the energy uncertainties cancel to first order and 
the effect of the losses in the forward beam pipe is diminished. Thus the optimal four 
"detector oriented" variables to characterize deep inelastic scattering at HERA are 
[EAS, 7 ]. 

Using these four input variables, there are three basic methods which make use of 
only two of them at a time, and which are precise enough to allow sensible kinematic 
reconstruction, namely: 

- The electron only method (e) which uses E and 9. 

- The hadrons only method (h) which uses £ and 7 [12|] 



- The double angle method (DA) which uses 9 and 7 JTO] . 

In the following, all methods (for instance the DA) using some information from the 
hadronic final state will be called "hadronic", although strictly speaking there is only 
one inclusive hadronic method (a complete set of formulae is given in the appendix). 
The h method will not be discussed in the following since it is not precise enough 
compared to the others, for neutral current DIS events. 

More than two variables are needed to determine the kinematics if the incident 
energy is unknown, which is the case when the incident electron emits a photon before 
the hard collision. This photon is often undetected since it is emitted colinearly with 
the incident electron beam direction, and can thus escapes inside the beam pipe. In this 
case three variables are needed to reconstruct the kinematics. For instance r/h = T,/2E a 
can be replaced by y^ = S/(S + S e ) to take into account the missing energy due 
to the escaping photon. A further important characteristic of y% is that at high y, 



1 The positive z axis is defined at HERA as the incident proton beam direction. 
2 We can define the similar quantities for the scattered electron: 
£ e = i? (1 — cos 9) p Tte = Esin6 i.e. tan | = ^. 
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E becomes the dominant term in the denominator and thus experimental errors on 
the E measurement tend to cancel between numerator and denominator. In the E 
method ]TTJ , Q| is constructed, like y% to be independent of QED initial state radiation 
(ISR) and to be optimal in terms of resolution; thus pT, e is used instead of PT,h'- 

Ql=P 2 T,e/^-y^) x s = Ql/sy s . (9) 

Insensitivity to ISR on x s is achieved simply by replacing s by 2 P D (E + E e ), thereby 
obtaining the IE method which is based on (E,6, E). The DA method was also rendered 
ISR independent by using E to reconstruct the initial electron beam energy. The IDA 
method [|lOj obtained in this way is thus based on (E,#, 7). The complete formulae 
and the comparison of these two methods can be found in [Tl||. They will not be 
considered in the following, since the gain obtained by the complete independence to 
ISR is not sufficient to compensate the loss of precision induced by the reconstruction 
of the incident electron energy. 

Fig. ^ shows the purity of the e, DA and E methods as a function of x in bins of 
Q 2 . The properties of the 3 standard methods to reconstruct the kinematics of neutral 
current DIS events at HERA are clearly visible: high precision of the e method at high 
y with a severe degradation at low y, good precision for the E method in the complete 
kinematic range, high precision of the DA method at medium (~ 10 2 GeV 2 ) and high 
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Figure 2: Purity at Q 2 =l, 5, 50, 500, 5000 and 25000 GeV 2 for the 3 kinematic 
reconstruction methods: e (open triangles), DA (open squares), E (closed circles). 
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Q 2 (~ 10 4 GeV 2 ) with a severe degradation at low Q 2 (~ 1 GeV 2 ). 



The "hadronic" methods have already been shown to display at low x and low Q 2 
a rather imprecise reconstruction of Q 2 . A simple solution to this problem is to use Q 2 
and to obtain, via Q 2 = xys, either x from a hadronic y or y from a hadronic x. For 
instance we have the mixed method |13[ (x m is obtained from y^ and Q 2 ) which has 
a good precision at low y IT5[. The precision at high y can be improved by using 



the mixed E (mS; y m s = y^ and Q^s = <3e)> or better the eE method (x e s = x^, and 



Qg S = Q e ) as was shown on fig. 2 and 3 of ref. [fll 



Another approach to combining two methods which have complementary properties 
has been tried by the HI collaboration in the analysis of the diffractive structure 
function |I1J. Since Q 2 = 4E 2 (1 — y)/tan 2 | both in the e and DA methods, an 
"average" method (labeled here ADA) has been introduced, in which yADA is obtained 
by weighting y e and yoA by y and (1 — y) respectively: 

= 2/ e + 2/£>a(1 - Vda) Q ADA = 7 ~e ' ( 10 ) 

tan | 

On the non- diffractive DIS events a better complementarity is actually achieved in the 
AE method by "averaging" the e and the E methods |T7]]. However neither the AS nor 
the ADA method bring an improvement compared to the simpler eE method as can 
be seen in fig. |3|: the AE method is giving similar performances to the eE one except 
at low x where it is slightly less precise, while the ADA method is better at high x and 
high Q 2 , but weaker elsewhere, in particular at low Q 2 . 
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Figure 3: Purity at Q 2 =l, 50 and 5000 GeV 2 for the 3 kinematic reconstruction 
methods: AE (open triangles), ADA (open squares), eE (closed circles). 

Before ordering these methods to underline their relations, let us introduce another 
method derived from the DA method, which will be useful later, for the understanding 
of more complicated methods and for the study of error propagation. Following the 
logic used for the E method, the inclusive hadronic angle 7 can be replaced by 7s 
defined by tan 7s/2 = E/pT,e since transverse momentum conservation implies pr,h = 



7 



PT, e even in the case of colinear ISR. This replacement improves the precision on the 
hadronic angle but has the obvious drawback of introducing a sensitivity to the electron 
energy reconstruction errors, which is absent in the DA method. This DE method 
defined using 7^ instead of 7 in the DA formulae, gives: 

Vds = tan7 S /2 = E/p T , e = ^ 
tan7 E /2 + tan. 0/2 E/p Ti£ + E e /p Tje 

The Q 2 DTj and xdt. are nevertheless different from their E counterparts (although xdt, = 
x m ) and we will see in the section 5 how the hadronic and electron miscalibrations affect 
them. 

In fig. f| the distribution of the reconstructed x (x rec ) divided by the true x (x gen ) 
and the equivalent distribution for Q 2 are compared for the e, DA, E, DE methods at 
high y (0.3-0.7)0. This figure, as all the figures in this paper, is obtained from a deep 
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Figure 4: Distribution of x rec /x gen (a) and Q 2 ec /Q 2 en (b) for the e, DA, E and DE 
methods for Q 2 > 7 GeV 2 at high y (0.3 < y < 0.7) 



inelastic sample simulated in detail in the HI detector (1995 set-up) as described in |] 
All known effects both from the detector and from the physics point of view (structure 
functions, QED radiation) are included in the simulation, which has been shown to 
give a good description of the HI data in the complete kinematic range [18[ |lp, |20 



The uncertainty on the electron energy, on the polar angle and on the hadronic energy 
scale are typically 1%, 1 mrad and 4% respectively. For figs. 4 and 5, a cut at Q 2 > 7 
GeV 2 has been made, which corresponds at HERA to the region measured with the 
highest precision. In fig. 4 the excellent resolution at low x and low Q 2 of the e method 
is clearly visible. Only the E method has a comparable, although lower, resolution in 
x. In Q 2 the weakness of all these hadronic methods is visible. The DE method is 
more precise than the DA one at low Q 2 , as expected from the way it is constructed. 
Unfortunately the e method is not precise at low y (as can be seen in fig. 2), so either the 



3 At low y, since all hadronic methods give rather similar results in terms of resolution, (contrary 
to the e method which strongly degrades), the comparison is better made using the plots of the purity 
(fig. 2). The purity of the DS method is almost identical to the S one at low y (not shown). 
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Table 1: Classification of the hadronic methods. For the methods using the Double- 
Angle approach, the correponding hadronic angle is also given. Note that the mS 
method is ISR independent only in y. The methods displayed in the last row are 
discussed in section 4. 

structure function measurement is done using different methods in different regions of 
the kinematic plane, or a precise method over the complete kinematic plane is found. 
Both these approaches have already been used at HERA, but for the more precise 
future measurements it is mandatory to focus on the second case and to optimize the 
reconstruction method to be used. 

To conclude this section, the hadronic methods are ordered in table [l|: each column 
represents the type of approach followed (Sigma, Double- Angle, Mixed), while each 
row represents a prescription to derive them: basic method, ISR independent, "opti- 
mized resolution". The relations become clear and the effects coming from the type 
of recontruction used can be disentangled, by comparing methods from the same row, 
from those due to the detector response which can be studied by comparing methods 
from the same column. 



4 Kinematic Reconstruction Improvements 

The description of the methods given above suggests that it might be possible to further 
optimize the reconstruction precision. The most direct way would be to perform a fit 
of the kinematics using all information available as described for instance in [f2~T| , |22f . 
However the redundancy is not large, since "p T h is not precise enough to be helpful in 
neutral current events, except at high Q 2 . Furthermore the method requires a good 
knowledge of the uncertainties of the measured variables over all the kinematic plane 
and the lack of statistics has prevented the production of a detailed enough map of all 
imperfections of the detectors. We will thus concentrate on the "analytic" improvement 
of the hadronic methods, in particular by "rescaling" the hadronic energy: S may be 
rescaled to approach its true value using the formula £ + £ e = 2 E D , if the error on E e 
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is assumed to be small compared to the E one (]. But if E is expressed as 2 E — E e in 
the h or E method, these methods become identical to the e method. An intermediate 
solution between the E and e method can be obtained by using the rescaling factor 
defined as 



2E n 



E + E e 

as was already implicitly done in some of the previously discussed variables: 



(12) 



ys = ry h ; y e x = r y s ; Ql = — ; Q 2 DTi =Q\r . (13) 

r 

An improvement on the E kinematic variables at low x is obtained by directly rescaling 
E, thereby defining an rE method: 

rE + E e 1 - y rS 

The corresponding "rescaled" hadronic angle is defined as tan7 rS = rE/p T e and allows 
to define a rescaled DE method (rDE). This method was already derived in a different 
way by the ZEUS collaboration, and is called the "PT" method^ |7j . Similarly an reE 
method is defined, using x r s and Ql. With these definitions we create the last row of 
table ^ with the "rescaled" methods. These new methods are compared to the e and 
E methods in x and Q 2 in fig. |5] in the same Q 2 ,y intervals as those of fig. |j. With 




the rescaling, xpt (avs) has indeed a better resolution than xdt, (^s)- However xpt 
is still less precise for y values between 0.3 and 0.7. than the simple, non-rescaled, xs, 
due to the propagation of the hadronic error in these two methods (cf section 5). The 

4 Indeed, the relative error on E increases when going to lower x and lower Q 2 since the hadronic 
final state contains an increasing fraction of low energy particles which are measured less precisely 
and which may go undetected in the HERA detectors. 

5 Actually, the "PT" method includes also a calibration of the hadronic final state using pr.e- Here 
we single out the analytic definition of the method. 
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difference in x between the rescaled methods and the e method is now smaller, which 
is also true in Q 2 where the rescaled methods (rE and PT) are similar and definitely 
better than the non-rescaled ones (S and DS). By using the Q 2 , the reE combination 
allows a further improvement with respect to the two other rescaled methods or to 
the eS method. The overall behaviour in the complete kinematic plane of one of the 
rescaled methods (the PT), compared to the eS can be judged by their purity shown 
in fig. ^| They both show a precise behaviour over the complete kinematic plane. The 
e£ method is more precise at low x, at high Q 2 , and has a comparable or better purity 
than the e method (compare fig. |2| and fig. even at low x and low Q 2 . 
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Figure 6: Purity at Q 2 =l, 5, 50, 500, 5000 and 25000 GeV 2 for the PT (open squares) 
and eS (closed circles) kinematic reconstruction methods. With the high precision 
measurement of the electron energy in the HI detector, the e£ method is better than 
the PT method at low x and also at high Q 2 (see fig. 7 for further comments on this 
matter). 

Before examining in the next section the propagation of the errors on the kinematic 
variables, let us consider the other characteristic of the rescaling, i.e. the increased 
dependence on the electron variables, since the rescaling factor depends also on the re- 
constructed energy and angle of the scattered electron. We can indirectly demonstrate 
the increased dependence by applying a further rescaling, and obtain three new meth- 
ods: r2S, r2DS and r2eS. The r2 rescaling factor is obtained recursively (assuming 
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tq — 1, r\ — r and y ro j] = ys) using the following formulae: 



and 



2 E 



(15) 



2/r„S = 



2/r n _iS + 1 - Ve 



r n _iS + S e 



The resolution at high y improves at each rescaling and actually y rn z — > y e when 
n — > ocQ implying that the three approaches (r n £,r n D£ and r n eS) converge to the e 
method. This means that the gain in precision at high y is obtained at the expense of a 
loss of precision at low y. Furthermore the "hadronic" method becomes more influenced 
by the e errors and the significance of the cross-check of the systematics which can be 
made when using the e and £ methods independently is reduced at each rescaling. 
However, we will now study the influence of the hadronic and electron reconstruction 
errors on these methods, to see when the usage of a rescaled method is worthwile. 

5 Error Propagation on the Kinematic Variables 

In this section we will first consider that the hadronic final state reconstruction is the 
dominant source of systematic error on the kinematic variables, and thus neglect the 
effect of the electron reconstruction, which will be considered specifically in the second 
part of the section. This approximation is legitimate at HERA, in particular at low 
Q 2 , since the energy carried by the hadronic final state which is visible in the detector 
becomes on average smaller as Q 2 decreases. Since the hadronic final state is a com- 
bination of several particles which enter the detector in different places, the precise 
effect of the detector response can be obtained only via a complete detector simula- 
tion/reconstruction program. However we can start by examining the propagation of 
the observable errors on the kinematic variables, before comparing the results to a 
realistic simulation. Namely we will consider errors like and SpT t h/PT,h which are 

related to the error on the hadronic angle £7 by [] 



The error on 7 arises from the spread of the hadronic final state particles in the detector, 
since it is not possible to ensure complete homogeneity in the detector response (even 
with a perfect absolute calibration, and with a perfect angular reconstruction, the 
fluctuations in the energy measurement have a different influence on £ and px,h as 
soon as there is more than one particle in the final state). In addition particle losses 
in the beam pipe affect £ and pT,h in a different way: losses in the forward (backward) 
beam pipe affect to a good approximation mainly pT,h (£)• This explains why, at low 
Q 2 , the DA method is not precise both at very low and very high y. The complete 

6 Since y r „ E can be written as 1+n+M 2 + .. + 1 tt „-i +tt „ /as , with u = 1 - y e 

7 In all the error equations, negative signs may appear in front of some terms, allowing to know the 
direction of the bias introduced by a given measurement error. When the terms are independent, the 
actual total error squared is obtained as usual, i.e. by quadratically summing the different terms. 



^7 



_ 6pr,h 
£ PT,h 



(16) 



sin 7 
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Table 2: Errors on y, Q 2 and x due to errors on the hadronic observables for the E, 
DS, eS, PT (=rDYl) and rS methods. For the coefficients in this table, the following 
abreviated notations are used: y = ; y r = y r s . 



error propagation of and ^7 for x, y and Q 2 for the h, DA and m methods are given 
in the appendix, but are not discussed further here since these methods do not allow 
a precise measurement over the full kinematic range. 

For the optimized and for the rescaled methods, the corresponding error propagation 
are given in table 0, noting that they depend only on one hadronic variable, S. At low 
y (y < 0.1), the errors of the three methods become identical: 

by 5E SQ 2 5x 5E 

y S Q z x S 

and a more precise study of the effect of the losses in the forward beam pipe and of 
the pT, e reconstruction uncertainties is needed since they become more important than 
those due to the hadronic energy uncertainty. 

At high y the error originating from the hadrons dominates and the situation is 
drastically changed: 

• For the optimized reconstruction methods, the S approach allows smaller errors than 
the DA or DS approach, as can be seen by comparing 5x^/x^ to Sxds/xdz in the y 
interval 0.3 — 0.7. The computed average ratio of these errors is 1/3 and is even smaller 
if comparing the S to the DA method. This explains the results already discussed of 

fig-!- 

• For the rescaled methods (PT and rE) the errors on x and Q 2 , with respect to 
the corresponding non- rescaled method are reduced by a factor (1 — y), providing an 
improvement at high y. Both methods have a reconstructed Q 2 which deviates from 
the true Q 2 by the same amount but in opposite directions. For x we have a similar 
situation as for the non-rescaled case since the improvement by a factor of (1 — y) 
in both methods does not cancel the difference in resolution observed above between 
is and xdt,, so the rS (and actually also the simple E) method provides a better 
estimation of x than the PT one. So, when the hadronic error is the dominant one, a 
method of the Sigma type should be chosen to have the best kinematic reconstruction 
over the widest kinematic range. 
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As already said, the "hadronic" methods described above are also influenced by the 
reconstruction of the electron angle and energy. Here we examine the situation at low 
Q 2 in which typical values for 5E/E and 56 in the HERA experiments are 1% and 1 
mrad. The case of high Q 2 events is treated in section 6. From the error propagation 
of E and 6 on the kinematic variables reconstructed with the e method, we can observe 
that the 6 error becomes significant only when 6 tends to its extreme (0 and 180°). 
At low Q 2 , i.e. if 6 is greater than about 175° its error dominates the total error. 
For the following discussion we do not consider 56 if it is not weighted by a tan 6/2 
factor large enough to render this error comparable to the error originating from the 
electron energy miscalibration. Considering only the uncertainty due to the angular 
measurement, we have: 

— ~0 — f — — tan- 00 — ~— tan- 00 18 

Ve Ql 2 x e 2 y ' 

For the DA method we have a different 6 dependence (see table 5): 

Uoa sin0 Qda srn ^ x da sin0 

implying a more severe degradation at low Q 2 than for the e method, since jjda is 
strongly affected by the error, and xda, Qda are more affected than x e , Q 2 . 

To study the other "hadronic" methods we observe that the errors on S e and r due 
to 56 are such that: 

5S e /S e ~ ; 5r/r ~ (20) 

implying that for all hadronic methods which use pT, e (S,DE, rS, PT, eS,..) the errors 
due to (see also table 6) are identical, namely: 

$y n 5Q 2 rn Sx 6 rn 

— ~0 -^r - - tan -56 — ~ -tan- 56 21 

y Q 2 2 x 2 y ! 

We thus have a simple situation: the error on propagates in the same way in the e 
and in all these hadronic methods and this fact comes from the usage of the transverse 
momentum of the electron. Since y is essentially not affected by the £0 error, similarly 
all hadronic y are also independent of this error. However, both Q 2 and x have a 
divergent error at large angle, implying that the measurement of the lowest Q 2 by 
using larger angles is hampered by the angle error propagation. This suggests the 
usage of lower incident (beam) electron energies to measure precisely the lower Q 2 
domain, rather than further increasing the acceptance at large angle as has already 
been achieved at HERA, by shifting the interaction position by about 70 cm in the 
forward direction and by upgrading the detectors close to the backward beam pipe. 

To study now the influence of the electron energy E on the standard and rescaled 
methods, we compare in table |3| the energy error propagation for the S, DS, eS, PT 
and rS methods, keeping in mind for comparison the well-known ^ behaviour 
of 5x e /x e . The lack of dependence on E in the DA method is changed to a weak 
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Table 3: Errors on y, Q 2 and x due to errors on the energy of the electron for the E, 
.DE, e£, PT and rE methods. The following notations are used: y = y-% ; y r = Z/rE- 



dependence in the DE method. The E dependence of Q 2 for the E method (and to a 
lesser extent for the e method: 5Q 2 /Q 2 = ^r), is stronger than the DE one. At high y 
however, the influence of E on the Q 2 of these three methods become similar. A similar 
effect is visible on x, for which the E dependence of the DE is small but increases at 
high y, while it is larger for the E and e methods at low y but decreases with increasing 
y. For the rescaled methods the same message holds, but the E dependence on the 
PT method is larger than on the DE one (a factor 2 — y larger, for all variables), and 
conversely, smaller for the r£ than for the E one (from oc2 — ytooc2 — 2y), thereby 
reducing the initial difference between the E and DE methods. 

As a general conclusion, the PT method is better suited for low y than for high 
y and is less influenced by the electron energy miscalibration. The usage of the PT 
method should thus be reserved to the case in which the miscalibration error on the 
electron is of the same order as the hadronic one. This can be verified by comparing 
fig. 0a, which shows the purity at Q 2 =20 GeV 2 of the e£ and of the three rescaled 
methods in the nominal calibration conditions of the HI experiment, to fig. [7|b, which 
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Figure 7: Purity at Q 2 = 20 GeV 2 for 4 kinematic reconstruction methods: r£ (open 
triangles), PT (open squares), eE (closed circles) and reT, (open circles), a) nomi- 
nal calibration; b) Electron energy smeared by 4%; c) Electron energy smeared and 
miscalibrated by 4%. 
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shows the same purities, but obtained after degradation of the electron energy response 
with a gaussian smearing of 4%, and to fig. 0c, which has an additional systematic shift 
on the electron energy scale of —4%. The figures show how the PT method remains 
relatively insensitive to the energy degradation, thereby becoming more precise than 
the other methods which are more affected by the degradation. 



6 Kinematics at High Q 2 

With the foreseen increase of luminosity, the HERA physics program in the next decade 
will focus on physics at large momentum transfer. These studies have been already 
started, and have been spurred by the observation of an excess of NC DIS events at 
high Q 2 (Q 2 > 15000 GeV 2 ) |23|, ||§. The difficulty of reconciling both the HI and 
ZEUS signals 0, ^] in the hypothesis of a single narrow resonance has underlined the 
importance of kinematic reconstruction at high Q 2 in particular when dealing with low 
statistics on which radiative effects can have a strong influence. For the future data, a 
simple comparison of the different methods can be made, based on the bin purities and 
using high simulated statistics. A more elaborate study would require for instance to 
study the optimal (x, Q 2 ) binning for each method. However the basic properties can 
be understood with a single binning choice well adapted to the goal pursued by the two 
collaborations, i.e. an uncertainty on the electron and hadronic energy measurement 
of about 1 and 2 % respectively. 

We compare in selected high Q 2 bins, in fig. || the optimal eS method to the two 
simple methods which are precise at high Q 2 , i.e. the e and DA ones. The DA method 
has a good purity in the whole (x, Q 2 ) range displayed, but is in general less precise 
than the e or eS method. These two last methods display a similar behaviour at high 
Q 2 , both at high and low x. At Q 2 below 10000 GeV 2 , the low y weakness of the e 
method starts to be visible in the highest x points (see also fig. 2). 
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Figure 8: Purity at high Q 2 for 4 kinematic reconstruction methods: e (open squares), 
DA (open triangles), eS (open circles) and u (closed circles). 
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These three methods are also compared to the on method which was introduced 
in II as an extension of the £ method in which the hadronic final state response is 
recalibrated using energy-momentum conservation, and assuming that the hadronic 
miscalibration effects propagate in a similar way in pT,h and S, i.e SpT t h/PT,h — 8H/H. 
This approximation is correct at Q 2 above a few thousands GeV 2 . For instance, the 
resolution of the inclusive invariant mass M u , which is related to x u via M= v / x • s, was 
shown to be excellent at high Q 2 . It can be seen in fig. |8], how this method is superior 
to the others at high x, while staying comparable to the best one at low x. Using Q 2 
instead of Q 2 U in the uj method (i.e. an "eu;" method, not shown) further improves its 
purity by about 5 % in all bins, rendering it optimal in all these high Q 2 bins. 

More generally, it can be noticed that the stability and purity are significantly 
higher at high Q 2 than at low Q 2 , thus a finer binning both in x and Q 2 could be 
used provided the recorded number of events is large enough. In this respect, the 
high luminosity HERA upgrade is indeed a very challenging program and will allow to 
explore in fine detail the quark-gluon structure of the nucleon at the highest Q 2 . These 
studies will also make use of the charged current events, for which the kinematics, 
which has not been studied here, also improves at high Q 2 (see for instance [Q). This 
is due to the improved precision on pr,h, itself coming from the stronger collimation 
and higher energy of the jet(s). 



7 Conclusion 



After a brief sketch of some experimental problems encountered in extracting the struc- 
ture function F2 from the measured cross-sections, we have reviewed and compared in 
detail the different ways of reconstructing at HERA the kinematic variables of the in- 
clusive deep inelastic scattering. At low x, the effect of rescaling the hadronic energy 
using both hadronic and electron observables has been studied and generalized to all 
methods, allowing one to derive new properties intermediate between the e and the 
S methods. The error propagations allow one to choose the most sensible method 
depending on the kinematic region under study and on the available precision of the 
energy and hadronic calibration. The PT method has been compared to the e, E, 
rS and eS method and is found to be well suited for low y but less precise than the 
other methods at high y (i.e. low x) if the electron energy reconstruction is precise 
(<5E/E ~ 1-2%). In this case the eS is the most precise method over the full kinematic 
range. At high Q 2 it has been shown that the kinematic reconstruction becomes very 
precise and that further optimization is obtained by using the complete information 
coming from the electron and from the hadronic system, as it is done for instance in 
the u method. This potential precision underlines the necessity to have high statistics 
to reach a systematic limitation on the measurement and renders the exploration of 
the dynamics of the high Q 2 collisions, which will be done with the HERA luminosity 
upgrade program, particularly exciting. 
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Appendix: 
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Table 4: y and Q 2 formulae for the electron and the main hadronic methods discussed 



in the text, using: E = E h (E h - p Z;h ), E e =E (1 - cos0), r = grpjr and r 2 = rS+ ^ 
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Table 5: Errors on y, Q 2 and x originating from the errors on the measured variables 
J2,PT,h, 7, E, and 9, for the h, DA and m methods. 
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Table 6: Errors on y, Q 2 and x due to the error on the polar angle of the electron for 
the e, S,DS and PT methods. In the eS method, y is two times larger than for the £ 

method, i.e. & = 2 ^4 50. In the rS method: % = - tan f + 2 ^ s(2 ~ ys) 50. 

' J/ tan £ Q 2 2 1 tan £ 
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